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Abstract 

The Dirac oscillator is an exactly soluble model recently introduced in the context 
of many particle models in relativistic quantum mechanics. The model has been 
also considered as an interaction term for modelling quark confinement in quantum 
chromodynamics. These considerations should be enough for demonstrating that the 
Dirac oscillator can be an excellent example in relativistic quantum mechanics. In this 
paper we offer a solution to the problem and discuss some of its properties. We also 
discuss a physical picture for the Dirac oscillator's non-standard interaction, showing 
how it arises on describing the behaviour of a neutral particle carrying an anomalous 
magnetic moment and moving inside an uniformly charged sphere. 

Resumen 

El oscilador de Dirac es un modelo exactamente resoluble que ha sido introducido 
recientemente en el contexto de la mecanica cuantica relativista de muchos cuerpos. 
El problema ha sido tambien explorado como posible fuente de un termino de inter- 
accion para modelar confinamiento en cromodinamica cuantica. Estas consideraciones 
establecen sin lugar a dudas que el oscilador de Dirac puede servir como ejemplo in- 
teresante en mecanica cuantica relativista. Este articulo ofrece una solucion al prob- 
lema y la discusion de algunas de sus propiedades. Tambien discutimos una imagen 
fisica que se ha introducido para el termino de interaccion del oscilador, mostrando 
como surge de considerar el comportamiento de una particula neutra con momento 
magnetico anomalo que se mueve dentro de una esfera cargada uniformemente. 



Classification numbers: ll.lO.Qr; 03.65.Ge; 12.40.Qq. 
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1. Introduction 

Interesting and exactly soluble problems in relativistic quantum mechanics 
(RQM) are scarse. Disregarding the one dimensional examples — of which the book 
of Greiner (1991) is probably the best source — the typical problems for which the 
Dirac equation is exactly soluble are (Schweber 1964) the hydrogen atom, a particle 
in an homogeneous magnetic field, a particle in the field of an electromagnetic 
plane wave, and the Morse oscillator. There seem to be no new additions to the 
traditional textbook set in the last few years. Though in some situations it would 
be desirable to discuss examples in RQM (not necessarily exactly soluble) more 
related with contemporary applications, most of these are too technical and this 
usually forbids the discussion. 

Our purpose in this paper is to discuss a problem in RQM, the so-called Dirac 
oscillator which, despite its various uses, has an straightforward exact solution. 
Though variants of the model were discussed some years ago (see, for example, Ito 
et al. 1967, Katriel and Adam 1969), the model was recently rediscovered in the 
context of relativistic many body theories (Moshinsky and Szczepaniak 1989). The 
Dirac oscillator has been also studied in connection with supersymmetric RQM, 
with quark confinement models in quantum chromodynamics (QCD) and with 
relativistic conformally invariant problems (Moreno and Zentella 1989, Benitez 
et al. 1990, Martinez- y- Romero et al. 1990, 1991, Martmez-y-Romero and Salas- 
Brito 1991). On the other hand, as we hope to make clear in this work, the problem 
has interesting features which can make it a helpful example in RQM. Besides the 
Dirac oscillator is closely related to the non-relativistic harmonic oscillator, as we 
exhibit in section 2. 

Our paper is organized as follows. In section 2 we introduce the Dirac oscil- 
lator and discuss its principal properties. In section 3 we offer a physical interpre- 
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tation for the potential of the Dirac oscillator, and cast its equation of motion in a 
manifestly covariant form. A complete solution for the problem is given in section 
4. We also show that the radial part of the components of its spinor eigenfunctions 
have the same form as the radial eigenfunctions of a 3-D non-relativistic harmonic 
oscillator. Conclusions and comments are given in section 5. 

2. The Dirac oscillator 

2.1. Relativistic quantum mechanics 

Before introducing the Dirac oscillator, let us first briefly recall the fundamen- 
tals of Dirac's RQM (Bjorken and Drell 1964, Greiner 1991). The basic equation 
in the theory is the free particle Dirac equation (for conciseness we use units such 
that h = c = 1): 

dip 

Hireeip = (a-p + mP) ijj = i— (1) 

where i^free is the free-particle Dirac Hamiltonian, p = — zV is the momentum, 
the 4x4 Dirac matrices are 

the a are 2x2 Pauli matrices and the I's and O's in P stand, respectively, for 
2x2 unit and zero matrices. If we introduce the gamma matrices 7° = /3, 7 = /3a, 
then, on multiplying by (3, the Dirac equation can be put in the manifestly Lorentz 
covariant form 

(i^^'d^-m)iP = 0. (3) 

An electromagnetic interaction is usually introduced in the free particle Dirac 
equation (1) (or (3)) using the standard minimal-coupling prescription p — > p 
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— eA, and H ^ H — e4>, or, in a Lorentz covariant fashion, p^^ ^ pi^ — eA'^, where 
A is the vector potential, cj) the scalar electromagnetic potential and = {(f), A). 
This prescription suffices for writing down the Dirac equation for most but not 
for all the problems, some are not amenable to this procedure. For example, for 
describing the electromagnetic interacion of a neutron, the Dirac equation must 
be written as (Bjorken and Drell 1964) 

{i^d^ + 4^^^-^'^' - ™n)0 = 0, (4) 

where F'^^ is the electromagnetic field tensor, e is the magnitude of the electron 
charge, k is the anomalous magnetic moment and is the mass of the neu- 
tron. As we discuss in the next section, the coupling term used for the Dirac 
oscillator can be interpreted as an anomalous magnetic interaction. Let us recall 
that anomalous interactions are simply a phenomenological way of describing the 
effect of residual electromagnetic interactions between electromagnetic fields and 
globally neutral composite particles with charged constituents. 

2.2. The Dirac oscillator Hamiltonian 

As a simplified model for complex interactions the harmonic oscillator has 
many uses in non-relativistic quantum mechanics, that it can be exactly solved 
is thus particularly useful. This fact is intimately related with its Hamiltonian 
being quadratic in both the momenta and the spatial coordinates. With these 
considerations in mind and trying to obtain an analogous model in RQM, the Dirac 
equation has been explored in the search of a potential that can make it linear 
in both the momenta and the spatial coordinates (Moshinsky and Szczepaniak 
1989); the resulting problem has been termed the Dirac oscillator. To obtain the 
non-trivial equation of motion for this oscillator, linear in both r and p, we have 
to perform the following non-minimal substitution in the free particle Hamiltonian 
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-^^f'ree 

p — > p — irmvPr, (5) 

where m is the mass of the particle described and co is the oscillator frequency. 
The Dirac equation for the system is then 

dip 

i— = Hip — (cK-(p — imiorP) + mP)ip. (6) 

At this point, we may wonder why the interaction is introduced as we have 
done and not by simply turning on a linearly growing potential in i^free? There 
are two reasons for this. First, the pure linear term can always be gauged away 
— to see this, introduce a linearly growing potential to H^ee and take A = as 
a gauge-function for recovering the free particle equation. Second, as we show in 
section 3, the inclusion of the /3 matrix is crucial for obtaining a Lorentz covariant 
interaction. We also remark that although the p — imujfiY term is obviously not 
hermitian, the complete Hamiltonian remains Hermitian due to the presence of 
the a matrix. The prescription (5) also guarantees the C, P, and T invariance 
properties of the Dirac oscillator (Moreno and Zentella 1989). In particular note 
that the parity invariance of iJ, i.e. the purely spatial invariance under r' = — r, is 
fairly obvious from equation (6) if you remember that in Dirac's RQM the parity 
operator is 

P = e*'^7^ (7) 

where e^"^ is an unobservable phase factor conventionally chosen as one of the four 
values ±1 or ±^ (Greiner 1991). 

2.3. Properties of the Dirac oscillator 
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We turn now to stablishing the relationship between the Dirac oscillator and 
the harmonic oscillator. Let us first take the square of the Hamiltonian of the 
Dirac oscillator, in this way, and after a few straightforward manipulations, we 
get 



H'^ = (a-(p - iujrp) + m(3f = p'^ + m^u?r^ + (4S-L - 3)r7j^/5, (8) 
where 



S = (9) 



is the spin and 



L = rxp (10) 
is the orbital angular momentum of the oscillating particle described. 

We can see now that, as equation (8) shows, becomes essentially a Klein- 
Gordon Hamiltonian with harmonic oscillator interaction plus a spin-orbit cou- 
pling term. In this sense, we may say that the Dirac oscillator is something like 
the "square root" of a linear harmonic oscillator. The squared Hamiltonian (8) can 
be used to obtain in a simple way the energy eigenvalues of the Dirac oscillator, as 
we show in section 4. We must remark that the squared Hamiltonian (8) is com- 
posed only of operators commuting with the /3-matrix — these operators are called 
even. This property of i?^ implies that a closed form of the Foldy-Wouthuysen 
(FW) transformation can be found for the problem. This useful point is discussed 
in Martmez-y-Romeroet al. (1990). 

If we define the total angular momentum of the Dirac oscillator in the usual 
way as J = L-|-o"/2, it is easy to show that the system conserves the total angular 
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momentum. To this end, let us show first that neither L nor a are separately- 
conserved: 



[L, B\ = i{axp) - m<jj{rXa)P, 



(11) 



and 



[<T/2,H] 



i{axp) + mu>{rxa)P, 



(12) 



where [A, B] = AB — BA. But, obviously, J = L + a/2 does commute with H 
and, therefore, the total angular momentum J is conserved by a Dirac oscillator 



3. Physical origin and Lorentz covariance of the interaction term. 

To obtain a physical model for the interaction term in the Hamiltonian of the 
Dirac oscillator we follow here the ideas of Moreno and Zentella (1989), see also 
Benftez et al. (1990). Let us begin with a simple electrostatic problem. Consider 
an uniformly charged dielectric sphere of radius i?, the electric field produced inside 
the sphere varies as E = — Ar (A is a constant), whereas the magnetic field vanishes 
everywhere B = 0. We can always consider a very large sphere to safely disregard 
any edge effects. We are going to show how a particle with an anomalous magnetic 
moment coupled to the electromagnetic field produced by the sphere leads to the 
Hamiltonian of the Dirac oscillator. In the process we also exhibit the Lorentz 
covariant properties of the latter. 

In the proper frame of the sphere (let us call this the laboratory frame), we 
may calculate the electromagnetic potential A^^ associated with its static electro- 



[J, H] = 0. 



(13) 
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magnetic field. Good expressions for this are given by A = 0, = Xtr, that 
is 

AU = mtr). (14) 

To put this expression in a Lorentz covariant form, we take advantage of the 
gauge invariance of the electromagnetic field. To this end we select the gauge 
function 

4^ 3 ^' 

to obtain the new potential 

<ab = <b-5'^A=^(^2 + r^2^r). (16) 

Introducing the unit four- vector Uj^^^) — (I5O) — which may be interpreted as the 
four-velocity of a point particle attached to the origin of the laboratory frame 
(Nunez- Yepez et al. 1989) — we can rewrite the 4-potential (16) in the form 

Ai'^^[2{U- x)x'' - x'^U"] . (17) 

This expression is manifestly Lorentz covariant, the electromagnetic field ten- 
sor produced by the sphere can now be calculated as 

F^, = XiU^x" - Wx^). (18) 

What is the relativistic quantum behaviour of a neutral quantum particle 
moving inside our dielectric sphere? In the laboratory frame only the electric field 
is non-vanishing and, at first sight, it would seem that a neutral particle could not 
interact with this field. But, as we already pointed out in section 2, this is not the 
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(15) 



case if the particle carries an anomalous magnetic moment. When we substitute 
our expression (18) for the electromagnetic field in the interaction term appearing 
in equation (5), we get 

= £A(i„.r), (19) 

to obtain an expression for the Hamiltonian in the Dirac equation, we only need 
to multiply equation (19) by the (3 matrix (Greiner 1991). What we get in this 
way is an interaction term which grows linearly with r. Hence, the Hamiltonian of 
this problem can also be obtained, as we did in section 2 for the Dirac oscillator, 
by means of a non-minimal coupling prescription: 

p ^ p - —enXrp. (20) 
2m 

We have to choose 

A=^, (21) 
en 

for reproducing the Dirac oscillator equation (6). That is, for obtaining a Dirac 
oscillator the charge density of the sphere, and hence the electric field, must be 
taken as proportional to the Dirac oscillator frequency u). We have thus exhibited 
that the Hamiltonian of the Dirac oscillator may be regarded as describing a 
neutral particle (a neutron for example) with anomalous magnetic moment and 
interacting with a static radially growing electric field. 

4. Eigensolutions for the problem. 

We now proceed to calculate the complete solution for the Dirac oscillator 
problem. We find that each component of the radial spinorial eigenfunctions have 
the same form as the radial wavefunctions of the 3-D non-relativistic oscillator. 
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Hence, the name oscillator applied to the problem may also be justified at this 
level. 

As we have seen in section 2, the total angular momentum J commutes with 
H and, therefore, the angular momentum is conserved in our system. From the 
rules of sum for angular momentum, we know that j = I ±1/2, where j is the 
total angular momentum and I is the orbital angular momentum quantum number, 
respectively. We also know that in our system parity is a good quantum number, 
so we use it to clasify the energy eigenfunctions. As the parity of the spinorial 
eigenfunctions of radially symmetric problems in Dirac's RQM is of the form (—1)' 
(Bjorken and Drell 1964), it is useful to define 



+1 if parity is (-)J+V2 
— 1 if parity is ( 



,-1/2. (22) 



in both cases, we have I — j + e/2. 



Since the Dirac oscillator conserves angular momentum, it is convenient to 
split the energy eigenfunction into a radial part and an angular part according to 
(± corresponds to the value of e) 



1 / nr)yfmiio.<p: 

ip{r,t) = -{ \eM-iEt), (23) 



r 



where yf^i and yj^i/ are the spinor spherical harmonics (Greiner 1991) 



(24) 



and the Yi^^ are standard spherical harmonics. Notice that due to the fact that 
the parity operator P (eq. 7) is proportional to /?(= 7°), yi has to be of opposite 
parity to yf, the only way to accomplish this is by taking I' = j — e/2. 
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To obtain the eigenfunctions and the energy spectrum of the Dirac oscillator 
we use (equation (8)). From that equation, we obtain that the differential 
equations governing the behaviour of both components of the wave function, are 
the sum of a harmonic oscillator equation plus a spin-orbit coupling term (this is 
another way to see that the Hamiltonian of the Dirac oscillator conserves angular 
momentum). Notice that for Z = j ' + e/2 we should have 

(j + l/2)(j + 1/2 + e) = l{l + 1), (25a) 
whereas, for Z' = j — e/2, we should have 

(j + l/2)(j + l/2-e)=r(r + l). (256) 

In this way, we have established that the square of the angular momentum satisfies 
(cf. the expression for (3 in equation (2)) 

L2 = (j + l/2)(j + l/2 + e/3). (26) 

From equation (26), we can get the contribution to the energy spectrum of 
the spin-orbit coupling term. Using the relation 

J2 + l2 + 2S-L+ ^, (27) 
and combining equations (26) and (27), we get 

S.L=-ie(2j + l)/3-i, (28) 

which only means that the spin-orbit coupling makes a constant {i.e. independent 
of r) contribution to the energy spectra. Besides, we know that for problems with 
spherical symmetry the square of the momentum is related to the square of the 
angular momentum by the relation (Greiner 1991) 
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valid for radial functions defined like in (23). From this equation and (26), we get 
for the square of the Hamiltonian of the Dirac oscillator 



E^-m =-—-^-\ \-m uj r -\- mij[e{2j -\- 1) - (30) 

This equation shows that the radial solutions of the Dirac oscillator must be 
the same as the solution for the non-relativistic 3-D harmonic oscillator (Benitez 
et al. 1990). This has to be so because, in (30), the term — equals the 
radial part of a non-relativistic harmonic oscillator equation plus the constant 
term mLo[e{2j + 1) - /3]. 

Let us now exhibit the cigcnfunctions of the problem. For the so-called "big" 
radial component of the wave function F{r), the energy is positive (Greiner 1991) 
and the orbital angular momentum must he I = j + e/2. Using the analogy with 
the harmonic oscillator, F(r) must be of the form 



Fn,i{r) = A{y/rriwry+^ exp(-ma;r^/2)iFi(-n, I + 3/2, mujr^); (31a) 

For the "small" component G{r), the energy is negative and the angular momen- 
tum is I' = j — e/2. Consequently, in this case the solution must be given by 



G'n,«'(r) = ^(V™^r)''+^exp(-rmjrV2)iFi(-n,Z' + 3/2,ma;r2). (316) 

In these equations, the iFi{a,b,c) are confluent hypergeometric functions 
(Davydov 1967, Gradshteyn and Ryshik 1980), A is a normalization constant and 
n = 0, 1, 2, . . . Employing the normalization 
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J ^t^dV = 1, 

the constant A may be easily evaluated, the result is 



(32) 



(n + / + l-i) +(n + /-|) 

(33) 

To obtain the energy spectrum we proceed as follows: we know that the 
spectrum of a 3-D harmonic oscillator is given by (Davydov 1967) 

EM = co(^N + ^y iV = 0,l,2... (34) 

where N = 2n + I is the principal quantum number. The energy spectrum of the 
Dirac oscillator is thus given by the contribution of the harmonic oscillator plus 
constant terms 

- = 2mEN + mu;[e{2j + 1) - P]. (35) 
Thus the energy spectrum is given by (Bemtez et al. 1990): 

E = {mu;[2{N + 1) + e(2j + 1)] + m^}^/^ (36a) 
for the positive energy states, and 

E = -{mu)[2{N + 2) + e(2j + 1)] + m^}^/^ (366) 

for the negative energy states. It is to be noted that every state with quantum 
numbers {N±s,j±s) for e — —1 have the same energy than a state with quantum 
numbers {N ± s,j ^ s) {s an integer). This is a consequence of the dynamical 
symmetry of the problem (Quesne and Moshinsky 1990). The spectrum is shown 
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n 



|2n+/-e/2+a/2 



(2n+2L + l-2e)!! 



in figure 1 where this property and a characteristic supersymmetryc pattern may 
be discerned (Bemtez et al. 1990). It is important to point out here the difference 
between this spectrum and the bound state spectrum of, say, a Dirac hydrogen 
atom. In the latter all the bounded levels are always described by admixtures 
of the big and the small components of the wave function, it would not possible 
then to employ for solving it the explicit separation we used for solving the Dirac 
oscillator problem (Berrondo and Mcintosh 1970, Benftez et al. 1990). In this 
sense, we may say that the Dirac oscillator is more akin to a free particle than to 
the Coulomb problem. 

As the above discussion was based in the square of the Hamiltonian of the 
Dirac oscillator and not in H itself, one may wonder if the results presented are 
correct or not. We can show that in fact they are indeed correct. We only outline 
the principal steps. To begin with, recall that for systems with spherical symmetry 
the a-p term of the Dirac equation — acting on functions defined like in (23) — 
can be written as 



where 



r V 2 



(37) 



I d 

CKj. = a-r; Pr = — T^f- (38) 

r or 

After substituting these expressions in the Hamiltonian of the Dirac oscillator 
we obtain, after some manipulations, the system of coupled differential equations 



d 1 

dr r 



e(j + -j) +mu}r 



G{r) = {E- m)F{r), (39a) 



and 
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— + - 

dr r 



F{r) = {E + m)G{r) 



(395) 



This system of equations can be solved using standard methods (Bem'tez et 
al. 1990), and the results are exactly the same as those given in Eqs. (27) and Eqs. 
(32). 



5. Conclusions. 

We have introduced and discussed the properties of a Dirac oscillator. We 
have shown that the Hamiltonian of the Dirac oscillator describes the interaction 
of a neutral particle with an electric field via an anomalous magnetic coupling and 
have found that the energy eigenfunctions are of the same form as those of the 
non-relativistic harmonic oscillator. The reason behind this is that the square of 
the Hamiltonian for the Dirac oscillator becomes essentially a harmonic oscillator 
plus constant terms. We have remarked that as the squared Hamiltonian of a 
Dirac oscillator is composed only of even operators then a closed form for the 
Foldy-Wouthuysen transformation can be found. 

We have mentioned also that this system has some remarkable symmetry 
properties. First, as it conserves angular momentum it must be 0(3) invariant, 
but in fact that it admits the larger dynamical symmetry algebra 5'0(3, 1)©S'0(4), 
producing a more degenerate spectrum than expected from 0(3) considerations 
only (Quesne and Moshinsky 1990). This can be seen in figure 1. The energy 
spectrum exhibits also a supersymmetric pattern (Benitez et al. 1990, Martmez- 
y- Romero and Salas-Brito 1991). This property and its relationship with the 
existence of an exact FW transformation is analysed in Martinez-y-Romero et 
al. (1990, 1991). The existence of a closed FW transformation implies that the 
positive and negative energy solutions never mix, independently of the intensity of 
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the interaction. This characteristic, shown by a class of systems besides this one 
(another example is a particle moving in an arbitrary, time independent, magnetic 
field), has been referred to as the stability of the Dirac sea ( Martinez- y- Romero 
et al. 1990). 

We hope to have made clear that the Dirac oscillator can be profitably used 
when modern examples are needed to illustrate the uses of the Dirac equation. 
It can be exactly solved with no more effort than a non-relativistic harmonic 
oscillator and many of its properties can be calculated and discussed without too 
much elaborate formalism, it can even be solved using purely algebraic techniques 
(Martfnez-y- Romero and Salas-Brito 1991, De Lange 1991a,b). Furthermore, its 
usefulness for modeling confinement in a Lorentz covariant way can make its study 
more motivating for some people. 
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Figure Caption 

Figure 1. 

The energy spectrum of a Dirac oscillator shown in two different ways. For con- 
venience, the numbers given correspond to {E'^ — inn?')/{'muj) instead of E. 

Notice that states with quantum numbers {N ± s, j ± s) for e = — 1 have the 
same energy than states with quantum numbers (iV±s, j =F'S) (s an integer). This 
is a consequence of the dynamical symmetry of the problem. 

There is also a typical supersymmetric pattern which is fairly clear in the 
squared spectra shown. The states with positive energy and e = +1 have the 
same energy than the states with negative energy and e = — 1 excepting for the 
constant term 2muj. 
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